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A MATHEMATICAL MODEL OF AN ACTIVE CONTROL LANDING GEAR 


FOR LOAD CONTROL DURING IMPACT AND ROLL-OUT 

John R. McGehee and Huey D. Carden 
Langley Research Center 

SUMMARY 

A mathematical model of an active control landing gear (ACOLAG) has been devel- 
oped and programed for operation on a digital computer. The mathematical model in- 
cludes theoretical subsonic aerodynamics; first-mode wing bending and torsional charac- 
teristics; oleo -pneumatic shock strut with fit and binding friction; closed-loop, series - 
hydraulic control; empirical tire force -deflection characteristics; antiskid braking; and 
sinusoidal or random runway roughness. The mathematical model was used to compute 
the loads and motions for a simulated vertical drop test and a simulated landing impact of 
a conventional (passive) main landing gear designed for a 2268-kg (5000-lbm) class air- 
plane. Computations were also made for a simply modified version of the passive gear 
including a series-hydraulic active control system. 

Comparison of computed results for the passive gear with experimental data shows 
that the active control landing gear analysis is valid for predicting the loads and motions 
of an airplane during a symmetrical landing. Computed results for the series -hydraulic 
active control in conjunction with the simply modified passive gear show that 20- to 30- 
percent reductions in wing force, relative to those occurring with the modified passive 
gear, can be obtained during the impact phase of the landing. These reductions in wing 
force could result in substantial increases in fatigue life of the structure. 

INTRODUCTION 

Dynamic loads in large aircraft resulting from landing impact and runway and taxi- 
way unevenness are recognized as a significant factor in causing fatigue damage and 
dynamic stressing of the airframe structure. In addition, ground-induced vibrations re- 
sult in crew and passenger discomfort and can contribute to reduction of the pilot’s capa- 
bility to control the aircraft. Such vibration problems have been encountered with some 
conventional transport aircraft and have required modification of the gear design to im- 
prove ride and handling qualities after the aircraft is in service. 

These ground-induced dynamic loads and vibration problems will be magnified for 
supersonic -cruise aircraft because of the increased structural flexibility of the slender- 


body, thin-wing designs and the higher takeoff and landing speeds. During design studies 
of supersonic transports conducted in the United States, investigations of the ground ride 
qualities of one particular design indicated extremely high vibration levels in the crew 
compartment during takeoff (ref. 1). The design philosophy for airplane oleo -pneumatic 
landing gears has been to obtain the lightest gear for the maximum designed sink rate. 

For supersonic-cruise transport aircraft, it may be necessary to concentrate on limiting 
the load applied to the structure by the gear to obtain satisfactory structural dynamic re- 
sponse characteristics and a satisfactory fatigue life. One potential method for improving 
ground operations of supersonic cruise aircraft is the application of active control tech- 
nology to limit the loads applied to the airframe by the landing gears. 

Analytical studies have been conducted to determine the feasibility of applying active 
controls to the landing gear during roll-out or taxiing over uneven runways or taxiways. 
The study reported in reference 2 indicated that a shock strut with a hydraulically con- 
trolled actuator in series with the passive elements of the strut (series -hydraulic control) 
possessed the most desirable dynamic properties and would be quite feasible to implement. 
Little published information is available for actively controlling loads applied to the air- 
frame by the landing gear during the impact phase. 

The purpose of this paper is to present and validate a mathematical model of an 
active control landing gear (ACOLAG) employing a series-hydraulic control for load 
control during impact and roll-out and to present results from the application of this con- 
trol concept to a simply modified main landing gear of a 2268 -kg (5000 -Ibm) class air- 
plane. The model is evaluated for predicting airplane loads and motions encountered 
during vertical-drop and landing-impact simulations for a rigid airframe with a passive 
(conventional) oleo-pneumatic shock-strut landing gear. The model is also used to deter- 
mine control parameters and airplane response for a series -hydraulic active control gear. 

Edwin L. Fasanella, of LTV Aerospace Corporation, reviewed the equations of 
motion and assisted in debugging the antiskid braking and active control subroutines of the 
computer program. 


SYMBOLS 

The units used for the physical quantities defined in this paper are given first in the 
International System of Units and parenthetically in the U.S. Customary Units. Measure- 
ments and calculations were made in the U.S. Customary Units. 

Aco area of control orifice, m^ (ff2) 

Aq area of opening in shock-strut orifice plate, m^ (ft^) 
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Ap cross-sectional area of metering pin, (ft^) 

A\ shock -strut hydraulic area (piston area), (ft^) 

A 2 shock -strut pneumatic area (cylinder area), m^ (ff2) 

JR aspect ratio of wing, b^/s^ 

a lift-curve slope at finite aspect ratio, (dCL,/da)^ 

aoo lift-curve slope for infinite JR, dC-j^ Jda 

b wing span, m (ft) 

^D,f aerodynamic friction -drag coefficient 

aerodynamic induced-drag coefficient 
^D,p aerodynamic profile -drag coefficient 

^d,co discharge coefficient for control servo valve orifice 

^d,o shock-strut orifice discharge coefficient 

Cl aerodynamic lift coefficient 

Cx numerical constant for each type of tire (this symbol is kj in ref. 3; tjrpe III, 

Cx = 0.53; type VII, Cx = 0.60) 

Cz vertical -force coefficient for each type of tire (type I, Cz = 0.02; types III 

and VII, Cz = 0.03) 

C |3 torsional damping coefficient of control sensor mass, N/rad-sec (Ibf/rad-sec) 
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C/3,cr 

Df 

Dt 

d 

^ce 

*^cf 

•^ct 

*^cw 

<^cwG 

deh 

dew 

dwwG 

dCL/da 
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critical torsional damping coefficient of control sensor mass, 

N/rad-sec (Ibf/rad-sec) 

fuselage aerodynamic drag force, N (Ibf) 
tail aerodynamic drag force, N (Ibf) 
wing aerodynamic drag force, N (Ibf) 
diameter of unloaded tire, m (ft) 

arm between airplane composite mass center and wing elastic axis in fuselage 
Xb-Zb plane, m (ft) 

arm between airplane composite mass center and fuselage mass center in 
fuselage Xb-Zb plane, m (ft) 

control -line diameter, m (ft) 

arm between airplane composite mass center and aerodynamic center of pres- 
sure of tail along Xb-axis, m (ft) 

arm between airplane composite mass center and wing mass center in fuselage 
Xb-Zb plane, m (ft) 

arm between airplane composite mass center and wing-gear interface in fuse- 
lage Xb-Zb plane, m (ft) 

arm between wing elastic axis and hub in body coordinate system, m (ft) 
arm between wing elastic axis and wing mass center along Xb-axis, m (ft) 
arm between wing mass center and wing -gear interface along Xb-axis, m (ft) 
lift-curve slope, rad”l 





arm between airplane composite mass center and wing-gear interface normal 
to Zb -axis, rn (ft) 

d2 distance between tire -runway interface and airplane composite mass center 

normal to Xg-axis, m (ft) 

d3 distance between airplane composite mass center and tail center of pressure 

parallel to fuselage longitudinal axis, m (ft) 

d4 distance between airplane composite mass center and fuselage mass center 

normal to fuselage longitudinal axis, m (ft) 

d5 distance between airplane composite mass center and wing center of pressure 

parallel to fuselage longitudinal axis, m (ft) 

d0 distance between airplane composite mass center and wing center of pressure 

normal to fuselage longitudinal axis, m (ft) 

E Young^s modulus for wing structural material, Pa (Ibf /ft^) 

e base of natural logarithms 

^A,ss total shock-strut axial force, N (Ibf) 

^G,X ground force at tire-runway interface acting along Xg-axis, N (Ibf) 

^G,Z ground force at tire-runway interface acting along Zg-axis, N (Ibf) 

^N,ss force normal to shock strut, N (Ibf) 

Fn nose -gear force along Zg-axis, N (Ibf) 

^ss, 1 shock -strut binding friction force, N (Ibf) 

^ss ,2 sliding friction force between shock-strut piston and cylinder, N (Ibf) 


5 



^T,N 

Ft,t 

^wB 

G 


g 

h 

la 

le 


Isect 


Iyy 

i 


KwB 

Kx,T 

K/3 


Lt 


L\v 


tire force normal to runway, N (Ibf) 
tire force tangential to runway, N (Ibf) 

wing bending force due to cantilever deflection of wing, N (Ibf) 

shearing modulus of elasticity for wing structural material. Pa (Ibf /ft 2) 
gravitational acceleration, m/sec^ (ft/sec^) 
height of wing box, m (ft) 

mass moment of inertia of wheel and tire about axle, kg-m^ (slug-ft^) 
mass moment of inertia of semispan of wing about elastic axis, 

kg -m2 (slug-ft2) 

area moment of inertia of wing chord section, m^ (ft^) 

pitching mass moment of inertia of one-half of airplane, with respect to plane 
of symmetry, about composite mass center, kg-m2 (slug-ft2) 

wing bending spring constant, N/m (Ibf /ft) 

wing torsional spring constant, N-m/rad (Ibf -ft/rad) 

tire fore-and-aft spring constant, N/m (Ibf/ft) 

spring constant between braking control sensor mass and wheel, N/m (Ibf /ft) 
aerodynamic lift force of tail, N (Ibf) 
aerod 3 Tiamic lift force of wing, N (Ibf) 
length of active control hydraulic line, m (ft) 
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Zgg shock -strut length, m (ft) 

length between upper and lower shock-strut bearings for fully extended strut, 
m (ft) 

1 2 length between lower shock-strut bearing and hub for fully extended strut, 

m (ft) 

ground-induced torque on wheel about axle, N-m (Ibf-ft) 
brake torque about axle, N-m (Ibf-ft) 

M^f torque about wing elastic axis applied to fuselage, N-m (Ibf-ft) 

Mew torque about wing elastic axis applied to wing, N-m (Ibf-ft) 

Mn nose -gear moment, N-m (Ibf-ft) 

m^ composite mass, kg (slugs) 

mf one-half of fuselage mass concentrated at fuselage center of gravity, 

kg (slugs) 

m^ mass assumed concentrated at hub (axle), kg (slugs) 

m^ semispan wing mass assumed concentrated at semispan wing c.g. on spanwise 

chord containing wing -geai' interface, kg (slugs) 

m^ braking control sensor mass, kg (slugs) 

N number of cycles to failure 

Nrc Reynolds number of airflow 

p tire inflation pressure. Pa (Ibf/ft^) 
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Pat 

Pec 

Ph 

Pec 

Po,a 

Pr 

Pi 

P2 

Rh 

r 

Sw 

s 

t 

Vac 

V2 

^eff 

^temp 


atmospheric pressure, Pa (Ibf/ft^) 

n 

effective active control pressure, Pa (Ibi/ft*^) 

hydraulic pressure. Pa (Ibf/ft^) 

control accumulator pressure, Pa (Ibf/ft^) 

absolute tire inflation pressure at zero vertical load. Pa (Ibf/ft^) 

rated inflation tire pressure. Pa (Ibf/ft^) 

hydraulic pressure in shock-strut piston. Pa (Ibf/ft^) 

Q 

pneumatic pressure in shock-strut cylinder. Pa (lbf/ft‘‘) 

amplitude of runway roughness, m (ft) 

deflected radius of tire, m (ft) 
planform area of wing, m^ (ft^) 

shock -strut stroke, m (ft) 
elapsed time after touchdown, sec 

volume of hydraulic fluid transferred by active control, m^ (ft^) 

volume enclosed by shock-strut cylinder, m^ (ft^) 

effective velocity of hydraulic fluid through control, m/sec (ft/sec) 

intermediate value of velocity of hydraulic fluid through control for computing 
pressure loss, m/sec (ft/sec) 
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Xb,Xg 

Xcp 

y 

Zb»Zg 

ZBe 

^er,b 

^ew,b 

^G,g 

zg 

Q! 

"e 

“i 

y 


maximum width of undeflected tire, m (ft) 

displacement along Xb~ or Xg-axis, m (ft) 

fore-and-aft shift of tire footprint center of pressure, m (ft) 

absolute value of spanwise distance from fuselage center line to chord con- 
taining wing -gear Interface, m (ft) 

displacement along Zb~ or Zg-axis, m (ft) 

bending deflection of wing elastic axis at spanwise chord of wing mass center 
relative to wing elastic axis at fuselage center line along Zb-axis, m (ft) 

relative displacement between wing elastic axis at spanwise chord containing 
wing mass center and wing elastic axis in fuselage Xb-Zb plane along Zb- 
axis, m (ft) 

displacement of wing elastic axis at spanwise chord containing wing mass cen- 
ter along Zb-axis, m (ft) 

ground elevation measured along Zg-axis, m (ft) 

initial height along Zg-axis between shock-strut fluid level and control- 
reservoir fluid level, m (ft) 

total angular displacement of wing chord at spanwise location of wing mass 
center relative to wing chord at fuselage center line, a\ + rad 

wing elastic rotation, rad 

angle of incidence of wing, rad 

angular displacement of brake sensor mass, rad 

ratio of specific heat of gas at constant pressure to that at constant volume 
specific weight of hydraulic fluid, N/m^ (Ibf/ft^) 


9 



A 


difference or change of value 


6 combined vertical tire deflection, m (ft) 

6q vertical tire deflection for pure vertical loading, m (ft) 

6y vertical sinking of tire due to fore-and-aft tire deflection, m (ft) 

^ slip ratio between wheel and surface 

?7 slope of runway surface, rad 

6 angle of pitch of airplane, rad or deg 

K tire pressure -rise parameter 

X-ce angle, in Xjj-Zb plane, between wing chord and line connecting wing elastic 

axis to airplane composite mass center, rad 

X-cf angle, in X^-Zb plane, between fuselage longitudinal axis and line connecting 

fuselage mass center to airplane composite mass center, rad 

X-cw angle, in Xb-Zb plane, between line connecting wing mass center to airplane 

composite mass center and line through composite mass center parallel to 
fuselage longitudinal axis, rad 

XcwG angle, in Xb-Zb plane, between line connecting wing -gear interface to airplane 

composite mass center and line through composite mass center parallel to 
fuselage longitudinal axis, rad 

Xgb angle, in Xb-Zb plane, between Zb-axis and line connecting wing elastic axis to 

hub, rad 

II friction coefficient between tire footprint and runway surface, or tire-runway 

friction coefficient (positive along negative Xg-axis) 

/ijj dynamic viscosity of hydraulic fluid, N-sec/m^ (Ibf-sec/ft^) 

jUgg I friction coefficient between shock-strut lower bearing and piston 
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MsSjU 

Pa. 

ph 

a 

tOQ 

^n,/3 

Subscripts: 

b 

c 

ce,b 

e 

f 

g 

h 

1 

max 

t-At 


friction coefficient between shock-strut upper bearing and cylinder 
mass density of air, kg/m^ (slugs /ft^) 
mass density of hydraulic fluid, kg/m^ (slugs /ft^) 
stress level, Pa (Ibf/ft^) 

angular displacement of wheel about axle, rad 

frequency of variation of sinusoidal runway elevation, rad/sec 

natural frequency of brake control sensor mass, rad/sec 


body -axis system 

airplane composite mass center 

between airplane composite mass center and wing elastic axis in body-axis 
system 

wing elastic axis 

fuselage mass center 

gravity-axis system 

hub mass center 

initial value 

maximum value 

value at previous time step 
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w wing mass center 

wG wing -gear interface 

X linear dimension along Xjj-axis 

z linear dimension along Z^-axis 

Dots over s3nmbols indicate differentiation with respect to time. 

MATHEMATICAL MODEL 

Variables considered significant in analyzing actively controlled landing gear sys- 
tems are shown in figure 1. Of these variables the following are included in this study; 
impact loads; runway-induced loads; braking loads; aerodynamic loads; flexible -airframe 
characteristics; and load isolation. Some examples of previously proposed landing gear 
active control concepts for load isolation are shown in figure 2. The series -pneumatic 
and the series -hydraulic controls have the actuator in series with the stiffness and damp- 
ing elements of the shock strut. The parallel -hydraulic control has the actuator in paral- 
lel with the stiffness and damping elements of the shock strut. In all these systems, con- 
trol is implemented by metering the working fluid (oil or air) to the actuator. For the 
present analysis the series -hydraulic control was selected on the basis of results of the 
feasibility study presented in reference 2. 

Three rigid-body degrees of freedom are included in the model: vertical transla- 
tion, fore-and-aft translation, and pitching about the airplane mass center. Since the 
model was devised to study active control concepts of a single main gear with pitch rota- 
tion included, lateral symmetry was assumed and a simplified representation (linear 
spring with no rebound) of a nose gear was included to absorb the rotational energy. Def- 
initions of axes and forces for the rigid-airframe option and for the rigid-fuselage — 
elastic-wing option are presented in figures 3 and 4, respectively. The model accounts 
for the following variables; sinusoidal or random runway unevenness; empirical tire 
force -deflection characteristics and antiskid braking; oleo-pneumatic strut with fit and 
binding friction; closed -loop, series -hydraulic control; first -mode wing bending and tor- 
sional structural -elastic characteristics; and theoretical subsonic aerodynamics without 
ground effect. These variables will be discussed in order in the following sections. 
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Runway Model 


Two general methods of analysis for runway unevenness have been considered in the 
study reported in reference 4: deterministic and statistical (power spectral). It was 
stated in reference 4 that the deterministic approach should be used for d3mamic analysis 
when the response characteristics of the airplane are nonlinear. Consequently, runway 
models included in the present analysis are a sinusoidal variation of runway elevation as 
a function of distance traversed along the runway or a random variation of runway eleva- 
tion as a function of runway distance. For the sinusoidal variation, the runway elevation 
is defined in terms of the frequency and the amplitude of undulation as follows: 

2G,gl|. = ^h + "G*) 


The slope of the runway surface at any instant in time is: 


^it 


= tan ^ 


•fG,g|t - ^G,g|t_At) 


Xh,gL At 


The random variation of runway elevation is defined as piecewise linear slopes in terms 
of runway distance traversed as follows: 

“G,g|, = 

For the random variation a table of runway elevation as a function of distance traversed 
is required as input. 


Tire and Brake Models 

The tire and brake models employed in this analysis are shown in sketch (a): 
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Empirical ti re force-deflection equation^.- The mechanical properties of the tire 
are defined from semiempirical equations given in reference 3. Equation (125) of refer- 
ence 3 is used to define the normal tire force during a landing impact as a function 

of tire deflection and physical properties of the undeformed tire. This equation is (in the 
present notation) 


Ft^N = 2.4 


^ - Cr 

w ^ 



where 


6 — 6q - 6y 


6v = 


0.01 


^T,T 

%,T 


%,T 


= Cxd(p + 4pj.) 



1/3 


The tangential tire force Fp^p during free rolling of the wheel is defined as that 
force required to decelerate the wheel such that the tire tangential velocity is equal to the 
hub velocity. The tangential tire force during skidding (slip ratios ? greater than zero) 
is defined as the tire -runway friction coefficient /lx multiplied by the normal tire force. 

Tire -runway friction coeff icient . - The tire -runway friction coefficients are modeled 
as a function of wheel-runway slip ratio ? (called in ref. 5) on the basis of experi- 
mental data presented in references 5 and 6. (See sketch (b).) The dashed line represents 
the model of the tire -runway friction coefficient at touchdown and during wheel spin-up. 

In the present analysis, the friction coefficient at touchdown is assumed to be 0.336 on the 
basis of data in reference 6. Eharing wheel spin-up the tire -runway friction coefficient is 
held constant until the wheel slip ratio decreases to a value of 0.06 (sketch (b)), which is 
the intersection of the braking and spin-up curves. The friction coefficient is then reduced 
along the braking -model curve to the value of zero friction for free rolling. Antiskid 
braking is initiated at the time of nose -gear touchdown, and the friction coefficient follows 
the braking curve to a slip ratio of approximately 0,25. The antiskid braking system then 
operates with a sensor mass acceleration corresponding to the slip ratio (approximately 
0.25) for maximum tire -runway friction coefficient. 
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Sketch (b) 


Braking model.- The antiskid braking system model presented in reference 5 was 
simplified in this paper by omitting the tire elasticity at the tire -wheel interface. The 
brakes are controlled by the angular acceleration of the control sensor mass, which is 
attached to the wheel through a spring -damper arrangement (sketch (a)). Brake applica- 
tion and release are controlled by preselected values of control-sensor -mass accelera- 
tions. The brake torque applied by the system is assumed to be torque limited; onset and 
offset torque rates may be independently changed, but vary linearly with time. 

Oleo -Pneumatic Shock -Strut Model 

A sketch of a typical oleo -pneumatic shock strut and a detailed description of the 
mechanics of the landing gear are presented in reference 7. During a landing impact the 
strut piston forces hydraulic fluid to flow through an orifice at a high velocity and creates 
a pressure drop across the orifice which resists the motion of the piston. (See fig. 2.) In 
addition to the hydraulic resistance to strut closure, the fluid flow into the strut cylinder 
increases the pneumatic (nitrogen) pressure, which also opposes strut closure. Internal 
bearing friction forces must also be considered along with the pressure forces. The total 
axial force developed by a passive (conventional) shock strut, with a metering pin, may be 
expressed as follows: 

Fa,ss = (Pi - P2)(Ai - Ap) + P2A2 + Tf7(Fss,l + Fgs,2)|tanh (2s)| 

|s| 
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where 


Fss,l - ^Mss,u + Mss,Z ^ ^)^N,ss 

This is a modified form of the equation presented in reference 7. Since the shock strut 
connects the airframe and the hub, the axial force imposes acceleration to the airframe 
while reacting to hub motion. 

Series -Hydraulic Active Control Gear Model 

The purpose of the series -hydraulic active control gear is to limit the force applied 
to the airframe at the wing -gear interface by regulating the damping force in the oleo- 
pneumatic strut. The model of the series -hydraulic active control gear is shown schemat- 
ically in sketch (c). The control consists of a servo valve, a high-pressure reservoir, a 
hydraulic pump, a low-pressure accumulator, an electronic control circuit, accelerome- 
ters, pressure transducers, and a linear potentiometer. 



Sketch (c) 


The following example illustrates the operation of the series -hydraulic active con- 
trol gear during a landing impact and roll-out. The mass at the wing-gear interface is 
assumed to remain constant during landing; thus, the acceleration of the wing -gear inter- 
face reflects the force applied at this position. For a specific airplane design and the 
measured sink rate, the limit force, defined from the limit acceleration, is determined as 
the product of the mass and the square of the sink rate divided by the available shock-strut 
stroke. A signal corresponding to the limit acceleration is input to the electronic control 
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circuit. The electronic control circuit (sketch (c)) continuously monitors the acceleration 
signal at the wing -gear interface and compares this signal with the limit acceleration sig- 
nal. The control is actuated when the wing -gear accelerometer signal exceeds the accel- 
eration limit signal by a preset tolerance. 

Figure 5 schematically illustrates three operational phases of the series -hydraulic 
active control gear during a landing impact. First, the gear contacts the runway at the 
measured sink rate and the shock-strut force builds as the strut stroking velocity s in- 
creases. The gear continues to operate in a passive mode (see fig. 5(a)) until the alge- 
braic sum of the forces at the wing -gear interface (airplane weight, aerodynamic forces, 
and shock-strut force) exceeds the control limit force. The second operational phase 
occurs during shock -strut compression when the electronic control signals the servo valve 
to move the spool to connect the piston to the low-pressure accumulator. The rate at 
which the spool is moved is determined within the electronic control as a function of the 
stroking velocity of the shock strut. Hydraulic fluid will then flow from the shock -strut 
piston into the low-pressure accumulator (see fig. 5(b)) to limit the force rise in the shock 
strut. As the shock -strut velocity decreases, the shock -strut hydraulic force will de- 
crease. If the wing -gear-interface acceleration drops below the limit acceleration and if 
the wing -gear -interface velocity is greater than a predetermined value, the rate of open- 
ing of the servo valve will be decreased to permit an increase of the hydraulic force and 
raise the shock-strut force. With continuing decrease in the shock-strut velocity and 
accompanying decrease in hydraulic force, the electronic control will command the servo 
valve spool to commence closing to maintain the wing -gear -interface acceleration within 
prescribed tolerances of the limit acceleration. When the servo valve spool returns to the 
fully closed position, the piston is isolated from the low-pressure accumulator. The third 
operational phase occurs as the wing-gear acceleration drops below the limit acceleration, 
and the electronic control commands the spool to connect the piston to the high-pressure 
hydraulic reservoir. Hydraulic fluid is then injected into the piston under high pressure 
to increase the hydraulic force. (See fig. 5(c).) The rate at which fluid is injected is 
determined by the rate of change of the shock-strut extension velocity. 

The transition from the iiupact phase to the ground-roll phase is accomplished by 
linearly decreasing the limit -force signal from the value used during the impact phase to 
a value to be used during the ground-roll phase. During this transition, the limit accel- 
eration tolerance will be increased to a value compatible with the active control gear de- 
sign. With the assumption that the mass of the airplane remains constant during landing, 
the momentum of the wing -gear interface is determined by the wing -gear -interface veloc- 
ity. The objective of the transition phase is to reset the signal for control limit accelera- 
tion and the tolerance signals to values compatible with a wing -gear -interface acceleration 
of zero for the ground-roll phase. This infers that the impact energy of the airplane has 
been dissipated. Transition to the ground-roll configuration is initiated when the wing- 
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gear -interface velocity, which is continuously monitored by the electronic control, be- 
comes equal to a signal input equivalent to an impulse divided by the airplane mass. For 
takeoff the control will be activated in the ground-roll mode. After takeoff, when the gear 
becomes fully extended, the control will be recycled into the impact phase of operation and 
will be deactivated during gear retraction. When the gear is extended for landing, the 
control will be activated in the impact mode. 

Control effects are mathematically simulated by modifying the passive-gear shock- 
strut pressure equations, which are 

Pneumatic 


P2 = P2,i 


^2,i 


.y 
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The principles used to modify the pressure equations are as follows: Ehiring shock-strut 
compression, the volume of incompressible fluid that normally would flow through the 
orifice in a small time interval is diverted to the low-pressure accumulator at a rate 
equivalent to the shock-strut compression rate and thus limits the hydraulic pressure in- 
crease. The modification of the passive -gear pressure equations for simulation of control 
effects is therefore based upon two control parameters: instantaneous volume flow rate 
into or from the piston and cumulative volume added to or removed from the piston. The 
pneumatic pressure is a function of the cumulative volume of fluid which has flowed from 
the piston to the control low-pressure accumulator or from the control high-pressure res- 
ervoir into the piston. Therefore, the equation for the pneumatic pressure incorporating 
active control is 


/ Y 

. A2S - V,ej ‘ 

The passive -gear shock-strut hydraulic pressure equation is modified to reflect control 
operation by solving for an effective shock-strut compression or extension velocity based 
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on the control flow rate, the area of the piston, and the instantaneous shock-strut velocity 
as follows: 


Pi = Ph 
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The control flow rate and incremental values of flow are determined through an 
iterative process. These values of flow rate and flow are used in the pressure equation 
to solve for new values of pressure and shock-strut force. The following procedure is 
used to determine the flow rate and flow during strut compression: The shock-strut force 

was controlled to limit the force at the wing -gear interface during the preceding interval. 
The state of the integrated variables has been determined for this shock-strut force. The 
shock -strut force is recomputed on the basis of the updated variables and prior shock- 
strut and control parameters. If the shock-strut force results in a force at the wing -gear 
interface greater than the control limit force, the control must lower the shock-strut 
force. The updated hydraulic pressure and the control accumulator pressure are used to 
compute a temporary velocity of flow through the control line and servo valve: 
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This temporary velocity is then employed to compute the effective pressure Pgg between 
the piston and the control accumulator by introducing line losses: 


Pec - Ph 


32|ijjvtgjjipZgf 

T2 

dcC 


rH(z£ 


s) 


The effective velocity of flow for the instantaneous shock-strut parameters is defined as 


Veff = - 


2(Pec ~ Pgc) 
PH 


1/2 


The area of the servo valve opening required on the basis of the change in shock-strut 
velocity over the computing interval is determined as 
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Aco t “ -^co t-At 


As (Aj - Ap) 
^djCo^eff 


The control flow rate is defined as 


^ac t “ ^d,co^co t^'^eff 


and the incremental flow occurring during the time interval is 


^ - Vac ^ 


The shock-strut force for the active gear is determined by the same equation used for the 
passive gear. If the shock-strut force again results in a wing force greater than the limit 
force, the control orifice area is increased and new values of flow rate and flow are com- 
puted on the basis of the shock-strut parameters at the beginning of the interval. The 
flow which occurred during the first iteration is removed from the cumulative flow and the 
new value of flow is algebraically summed. This iterative process continues until the 
wii^ force is within the design tolerance of the control limit force. Mathematical simu- 
lation of flow from the control high-pressure reservoir into the shock-strut piston is 
accomplished in a similar manner. 


Elastic -Wing Model 

Limited structural -flexibility effects are considered by including in the mathemati- 
cal model first-mode wing bending and wing torsional characteristics. The wing bending 
mode is defined as the deflection of the elastic axis of the wing deforming as a cantilever 
beam. The wing mass, wing -gear interface, and engine mass are assumed to be located 
on the wing chord at the same spanwise location. The wing is also assumed to be rec- 
tangular in planform. The equation for the spanwise bending spring constant may, with 
the foregoing assumptions, be expressed as 

^ _ SETsect 

^wB o — 

The bending deflection zgg of the wing is determined from the relative motion of the in- 
tersection of the wing elastic axis with the plane containing the fuselage center line 

and the wing elastic axis at the spanwise chord containing the wing mass center. The wing 
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bending force is determined from the product of wing bending deflection and the wing bend- 
ing spring constant: 


FwB “ J^wB^Be 


The first-mode torsional characteristics of the wing are determined assuming a 
constant wing cross section from the fuselage center line to the spanwise chord containing 
the wing mass center and the wing -gear interface. Wing torsional deflection is assumed 
to occur about the elastic axis and the wing chord is assumed to remain undeformed. With 
these assumptions the wing torsional spring constant is 


^WT 
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where 

Ssect chordwise area of wing cross section, m^ (ft^) 

y spanwise distance from fuselage center line to chord containing wing-gear 

interface, m (ft) 

h height of wing box, m (ft) 

twbh thickness of horizontal wing box units, m (ft) 

twbv thickness of vertical wing box units, m (ft) 

Wwbx width of wing box, m (ft) 

The torsional elastic deformation ofe of the wing chord at the spanwise location of the 
wing mass center relative to the wing chord at the fuselage center line is determined from 
the relative angular motion of these sections. The wing toi'que is defined as the product 
of the torsional elastic deformation and the wing torsional spring constant: 

^ew “ ^wT^e 
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Aerodynamic Model 


The aerodynamic model employed is based on theoretical subsonic aerodynamics, 
neglecting ground effects, determined from reference 8, The theoretical lift -curve slope 
for a wing of infinite aspect ratio is 

aoo = dCL /da = 27 t per radian 


For a finite aspect ratio, the lift -curve slope is approximately 


a = 


aoo 



per radian 


The wing lift force is 

= O,5PaSwX0^gCj_, 


where 


Cl = a(0 + O' + tan"^ 

\ ^c,g 

The tail lift force is found in a similar manner. 

The wing aerodynamic drag force (ref. 8) is assumed to consist of induced drag, 
friction drag, and profile drag and is defined as follows: 

D\v ~ 0‘5paSw^c,g(CD,i + ^D,f + ^D,p) 


where 

Induced drag, 
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Friction drag, 

Co,f = 0.0375(NRe)'°-^^ 
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Profile drag, 


CD,p = 0.075(NRe)-°-^^ 

The drag forces for the tail surfaces and the fuselage are determined in a similar manner. 

EQUATIONS OF MOTION 

The deterministic solution consists of determining the forces to be used in the equa- 
tions of motion, defining the coupled equations of motion which govern the response of the 
airplane, and performing numerical integration of the equations of motion on a digital com- 
puter. The equations of motion were written for two separate options: rigid-airframe 

option, and rigid-fuselage — elastic-wing option. The coupled equations of motion are pre- 
sented in appendix A for both options. 

For the elastic wing, the wing mass center has additional degrees of freedom, which 
are perturbations in the Z^-direction and in rotation about the elastic axis relative to the 
rigid-body motions of the airplane. 

The equations of motion were solved by performing a numerical integration. The 
numerical integration procedure is based on the classical fourth-order Runge-Kutta form- 
ula, and use of Richardson's extrapolation-to-the -limit theory results in a fifth-order 
integration procedure. This procedure operates with a variable time interval which is 
sized to meet a specified local relative truncation error. 

ANALYTICAL MODEL SIMULATIONS 

To validate the mathematical model described in the foregoing sections and to per- 
form an evaluation of the active -control gear, a landing-gear and airplane configuration 
was modeled to permit simulations of vertical-drop tests, landing impacts, and landing 
roll-outs with the ACOLAG computer program. To accommodate the active control, the 
passive gear was modified in the analysis to provide additional stroke capability. The 
stroke was increased from 0.187 m (7.375 in.) to 0.264 m (10.375 in.) by increasing the 
cylinder length, pneumatic volume, and piston length. 

Passive Gear 

Vertical -drop test . - To permit comparison between computed data from the ACOLAG 
computer program with experimental drop -test data, the program was adapted to include 
tire hysteresis for the nonrotating wheel (ref. 3), and braking simulation was excluded. A 
lumped mass at the gear trunnion was used to simulate a rigid airframe. The passive 
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landing gear (ref. 7) from a 2268-kg (5000-lbm) class airplane was selected, since exper- 
imental vertical -drop -test data were available for this gear. A sketch of the shock strut 
and the geometric characteristics of the gear are presented in appendix B. The impact 
parameters were a sink rate of 2.7 m/sec (8.8 ft/sec); a lift force of 10.7 kN (2412.5 Ibf), 
which was held constant throughout the impact; a pitch ar^le of 0° (gear normal to impact 
surface); and a wheel rotational velocity of 0 rad/sec (0 ft/sec). 

Landing impact. - For comparison of computed data with data obtained from trans- 
versely symmetrical airplane landings, tire hysteresis effects were eliminated (ref. 3) 
and antiskid braking was included. An airplane configuration was defined in terms of 
geometry, a rigid -airframe lumped mass distribution, and aerodynamics compatible with 
the selected landing gear. The geometry and lumped mass distribution of the airplane 
configuration are shown schematically in figure 6. 

Since experimental landing -impact data are not available for the selected passive 
gear (ref. 7), ACOLAG was evaluated for landing impact by comparing the results of this 
simulation with computed results obtained from a multiple -degree -of -freedom takeoff and 
landing analysis (TOLA), which, as described in reference 9, will only accommodate con- 
ventional gears and a rigid airframe. The impact parameters were a sink rate of 2.7 
m/sec (8.8 ft/sec); a ground speed of 45.7 m/sec (i50 ft/sec); a pitch angle of 10.5°; an 
angle of attack of 13. 8°, which results in a lift force equal to the airplane weight at touch- 
down; and a wheel rotational velocity of 0 rad/sec (0 ft/sec). 

Landing roll -out. - To determine the damping characteristics of the modified passive 
gear, a landing roll-out simulation was made. For this simulation the airplane was placed 
on the runway with a pitch angle of 0 rad (0 deg), the gear and tire deflected to support the 
mass of the airplane, the tire and wheel in the free-rolling condition, and a ground speed 
equal to 44.5 m/sec (146 ft/sec). The simulation was started with the tire encountering a 
sinusoidally varying (3 Hz with a half -amplitude of 0.0254 m (0.0833 ft)) runway bump 
which had an initial elevation of zero. At the end of one -quarter cycle, the surface was 
dropped vertically to zero elevation and remained at this elevation for the remainder of 
the simulated roll-out. 


Active Gear 

Vertical -drop ^st. - No experimental data are available for the active control gear. 
Therefore, evaluation of the active control gear was accomplished by comparing computed 
results for the active gear with those of the passive gear for the vertical -drop test. The 
gear parameters and impact parameters were the same as those for the passive gear. 

Landing impact. - For this case evaluation of the active control gear was accom- 
plished by comparing computed results for the active gear with those of the modified pas- 
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sive gear. The airplane configuration and impact parameters were the same as those de- 
scribed in the section ” Landing impact’’ for the passive gear. 

Land ing roll-out . - To compare the damping characteristics of the active gear with 
those of the modified passive gear, a landing roll-out simulation was made for the airplane 
and runway parameters cited in the section ’’Landing roll-out” for the passive gear. 

The modification of the passive gear discussed in this paper is simply an expedient 
method to characterize the range of control parameters that would be required for such a 
design. However, it should be noted that an active control gear would not be a simple 
modification of a passive -gear design, but would require an independent design based upon 
optimization of gear parameters in conjunction with control parameters. 

RESULTS AND DISCUSSION 

To evaluate the ACOLAG computer program, results were obtained from the 
vertical-drop -test simulation and compared with the experimental data reported in refer- 
ence 7. The TOLA computer program (ref. 9) was also employed to simulate the same 
vertical-drop test and the results were compared with the experimental data. Evaluation 
of the ACOLAG landing -impact simulation was accomplished by comparing results obtained 
by using the rigid -airframe and passive-gear modes with results obtained for a similar 
calculation employing the TOLA computer program. 

Results are presented for vertical-drop-test simulations with stationary and verti- 
cally oscillating impact surfaces; for landing-impact simulations on a flat runway surface 
and a sinusoidally varying amplitude surface; and for a landing roll-out on a runway sur- 
face having a single, quarter -cycle, sinusoidal bump ending in a vertical discontinuity at 
peak amplitude, followed by a flat surface. 

Passive Gear 

Vertical -drop simulations. - Figure 7(a) presents computed and experimental time 
histories of shock-strut force and stroke during the impact phase of the vertical-drop 
test of the passive gear. Since rebound of the upper mass (airframe mass) occurred at 
approximately 0.16 sec in both the experiment (ref. 7) and the computer simulation, data 
are not presented beyond this time. The shock-strut forces computed from ACOLAG and 
TOLA are in excellent agreement and both are in good agreement with the experimental 
data. As was the case with the shock -strut force, the calculated shock -strut strokes for 
the two programs are in excellent agreement and both are in good agreement with the ex- 
perimentally determined strokes for most of the time history. 

Time histories of vertical ground force and tire deflection, computed for a rigid air- 
frame and a passive gear, are compared in figure 7(b) with experimental data obtained 
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from the vertical -drop test (ref. 7), Computed vertical ground forces are generally in 
excellent agreement with the ground force obtained from the experiment. The ground 
force computed from the ACOLAG program is in excellent agreement with that computed 
from the TOLA program. The tire deflections calculated from the ACOLAG program are 
in excellent agreement with those computed from the TOLA program and both are in good 
agreement with experimental data. 

Figure 7(c) presents a comparison of mass-center accelerations (nomenclature for 
experimental data is upper -mass accelerations) for the computed results and experimen- 
tal data. Since the shock-strut force is the only varying force applied to the upper mass, 
the trend of the comparison should be the same as that for the shock-strut force (fig. 7(a)). 
The trend is the same and the agreement between computed and experimental mass -center 
accelerations is good. 

Landing-impact simulations . - Evaluation of the ACOLAG computer program for the 
landing -impact option is accomplished by comparing computed results from the ACOLAG 
program with results from the TOLA program. These comparisons are presented in 
figure 8. 

Time histories of pitch attitude are compared in figure 8(a) for a landing impact on 
a flat runway. The agreement between the computed results from the two programs is 
excellent. 

Time histories of shock -strut force and stroke calculated from the two programs 
are compared in figure 8(b). The shock-strut forces are in excellent agreement. Al- 
though the shock-strut strokes computed from the ACOLAG program are slightly less than 
those calculated from the TOLA program, the agreement is good. 

Computed time histories of vertical ground force and tire deflection are compared 
in figure 8(c). The vertical ground forces and tire deflections calculated with the ACOLAG 
program are in good agreement with those computed with the TOLA program. 

Airplane mass -center accelerations obtained from the ACOLAG and TOLA programs 
are compared in figure 8(d). The Zb -axis accelerations are in good agreement to a time 
of approximately 0.17 sec. At this time, the TOLA program indicates that the nose gear 
has contacted the surface, and hence the Zb-axis accelerations increase. The ACOLAG 
program has the nose gear contacting the surface at a time of approximately 0.18 sec, at 
which time, because of the simplified nose -gear representation, data from the ACOLAG 
program are not comparable with those of the TOLA program. Mass-center Xb-axis 
accelerations are in good agreement to a time of approximately 0.11 sec. At this time the 
wheel has spun up and is in a free -rolling state. When the wheel is freely rolling, the 
tire-runway friction coefficient in ACOLAG is decreased to a value which produces no 
wheel slippage. This tire -runway friction coefficient resulted in a very small component 
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of ground drag force applied at the mass center, while the X^-axis component of the ver- 
tical ground force was sufficiently large to develop a positive X]^-axis acceleration of the 
mass center. The definition of the tire-runway friction coefficient in TOLA resulted in a 
larger tire -runway friction coefficient during free rolling of the wheel. Consequently, the 
ground drag force applied at the mass center and hence the mass -center X^^-axis acceler- 
ations were greater. Comparison of Xj^-axis accelerations was also terminated at the 
time of contact between the nose gear and the surface. 

Considering the good agreement shown between TOLA and ACOLAG results when 
compared with experimental data and when compared with each other, the ACOLAG pro- 
gram appears to be a valid tool for the study of various actively controlled landing gears. 

Series -Hydraulic Active Control Landing Gear 

Results which demonstrate the operation of a series -hydraulic active control gear 
for vertical -drop and landing -impact simulations are presented and compared in figures 9 
to 13. 

Vertical -drop simulations. - Computed results for the passive and active control 
gears are compared in figure 9 for a vertical drop onto a flat stationary surface. Fig- 
ure 9(a) presents a comparison between the wing forces experienced during a vertical- 
drop simulation of the passive and active gears. The wing force is shown as a function of 
wing displacement. Presentation of the active -gear results is discontinued at the maxi- 
mum value of wing displacement because in the vertical -drop simulation the aerodynamic 
lift was maintained constant at a value equivalent to the mass of the airplane, and hence 
results calculated during rebound are not physically realistic for a landing simulation. 

For the active-gear simulation the wing force was limited to 21.46 kN ± 334 N (4825 Ibf 
± 75 Ibf), which represents a 26 -percent reduction from that of the passive-gear 
simulation. 

It was noted in reference 2 that for aluminum structures undergoing a cyclical level 
resulting in a fatigue life of approximately iO^ cycles, the relation between the number of 
cycles to failure N and peak-to-peak stress a is 

N =: A/aS 

where A = 1.8 x lO^'^ Pa^ If a fully reversed stress of 241 MN/m^ (35 000 

psi) for the passive gear simulation is assumed on an aluminum wing structure, the fatigue 
life of this structure would be 34 000 cycles. With the wing force controlled to a value of 
74 percent of the passive -gear value, the same structure would have an increased fatigue 
life of approximately 4.5 times, or 154 000 cycles. Thus it appears that small reductions 
in cyclical stress levels have a significant effect on the fatigue life of the structure. Re- 
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duction in wing force, however, results in an increase in wing displacement of approxi- 
mately 10 percent for this case, which is primarily reflected in an increased shock-strut 
stroke. For energy compatibility to exist between the passive -gear and the active -gear 
simulations, the hatched areas of figure 9(a) must be equal. The energies represented by 
the hatched areas were determined to be nearly the same. Thus, the incorporation of the 
control analysis into the passive -gear analysis results in an energy -compatible procedure, 
which increases confidence in the validity of active control landing gear analysis. 

Time histories of shock-strut force and stroke for the passive and active gears are 
compared in figure 9(b). At approximately 0.047 sec, control of the shock-strut force 
level at 21.46 kN ± 334 N (4825 Ibf ± 75 Ibf) was initiated to maintain the wing force 
within the preset limits, which were exceeded at this time. Since the strut was compress- 
ing, the control removed fluid and the maximum flow rate attained was -428 £/min (-113 
gal/min) at 0.084 sec. A transfer from the low-pressure mode to the high-pressure mode 
was required at 0.11 sec to maintain the wing force. This transfer results in a shift to the 
lower bound of the control limit force and accounts for the lower shock-strut force during 
the time from 0.11 sec to 0.143 sec. The maximum injection flow rate of 458 f/min (121 
gal/min) occurred at a time of 0.129 sec. The momentum of the upper mass becomes 
equal to the impulse designed into the electronic control circuit at a time of 0.143 sec, and 
the electronic control commands a linear decrease in the limit force from the impact limit 
force to the roll-out limit force. When the roll-out limit force is reached, data presenta- 
tion is terminated. At the time of initiation of control at the roll-out limit force (0.18 sec), 
the shock-strut stroke for the active-gear simulation was approximately 17 percent greater 
than the stroke shown for the passive -gear simulation. 

Figure 9(c) shows calculated time histories of vertical ground force and tire deflec- 
tion for the vertical-drop simulations of the passive and active gears. As would be ex- 
pected for the controlled lower shock-strut force, both the vertical ground force and tire 
deflection were lower for the active -gear simulation than those for the passive -gear 
simulation. 

Computed results for the passive and active control gears are compared in figure 10 
for a vertical drop onto a vertically oscillating surface. The vertically oscillating surface 
was introduced to simulate a sinusoidally varying runway -unevenness effect, since no 
ground speed is involved with vertical -drop testing. The surface oscillation, a frequency 
of 10 Hz and a half -amplitude of 0.0254 m (0.0833 ft), was selected as a worse-case con- 
dition, since this corresponded to the resonant frequency of the hub mass. At touchdown 
the surface was set into an upward motion opposite to that of the gear. 

Wing force as a function of wing displacement is shown in figure 10(a) for the 
vertical-drop simulation. In the initial phase of this simulation the ground motion passed 
through the peak displacement and was moving away from the tire. The resulting initial 
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rise in wing force did not exceed the upper bound of the limit force (21.80 kN (4900 Ibf)); 
consequently, the control was not activated. During the fourth quarter of the first cycle 
of ground motion, however, the control was activated because the ground was again mov- 
ing into the tire and the resulting wing force exceeded the upper bound of the control limit 
force. The wing force, for the active gear in this simulation, attained a value of only 64 
percent of that of the passive gear with an increase in wing displacement of only 5 percent. 

Time histories of shock -strut force and stroke are shown in figure 10(b) for the 
vertical-drop simulation with the oscillating surface. The appreciable reduction in wing 
force (36 percent) with the relatively small increase in wing displacement (5 percent), as 
shown in figure 10(a), was not attained without penalty. To achieve these values, the max- 
imum control flow rates required, as shown in figure 10(b), were -723 ^min (-191 
gal/min) during compression and 1423 H/min (376 gal/min) during transition from impact 
to roll-out control. Also, a 25-percent increase in shock-strut stroke above that of the 
passive gear was required. However, it should be recalled that the surface motion was 
chosen to create a worse -case condition, and for this study the gear was not optimized to 
be compatible with the series -hydraulic control. 

Figure 10(c) shows time histories of vertical ground force and tire deflection com- 
puted for this simulation. As would be expected, the maximum ground force and tire de- 
flection for the active gear were smaller than those computed for the passive gear during 
the portion of the impact phase prior to rebound. 

Landing -impact simulations. - Comparisons of computed data for the modified pas- 
sive gear and the series-hydraulic active control gear are presented in figures 11 and 12. 
Figure 11 shows comparisons of computed data for landing -impact simulations onto a 
smooth runway. Comparisons of computed data for landing -impact simulations onto a 
sinusoidally uneven runway (frequency of 10 Hz with a half -amplitude of 0.0254 m (0.0833 
ft)) are shown in figure 13. 

Figure 11(a) presents comparisons of computed wing force as a function of wing dis- 
placement for the passive gear and the active gear. The active gear achieved a 20 -percent 
reduction in wing force relative to the passive gear, with an increase in wing displacement 
of 8 percent. 

Figure 11(b) compares computed time histories of shock-strut force and stroke for 
the passive and active gears. The shock -strut force for the active control gear gradually 
increases during the initial control phase, in contrast to the essentially constant shock- 
strut force for the vertical-drop simulation. (See fig. 9(b).) Since the wing force is the 
force being maintained by the control, the shock-strut force increases to offset the de- 
crease in lift force as the airplane pitches down during the impact. The maximum control 
flow rates required during this simulation were 447 C/min (118 gal/min) for injection of 
fluid into the piston and -397 ^/'min (-105 gal/min) for removal of fluid from the piston. 
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At the time of nose gear contact (0.18 sec), the shock -strut stroke of the active gear was 
25 percent greater than that of the passive gear. 

Time histories of vertical ground force and tire deflection are shown in figure 11(c). 
Again, the maximum vertical ground force and tire deflection for the active gear were less 
than those of the passive gear during the impact phase. 

Figure 12 presents time histories of control flow rate and control flow which were 
required during the landing -impact simulation on a smooth runway. Following control 
initiation, the control removed fluid from the strut and attained a maximum flow rate of 
approximately -360 £/min (-95 gal/min). The flow rate decreased from this maximum 
value to zero at a time of 0.098 sec, which indicates that the control required a transition 
from the low-pressure mode of operation to the high-pressure mode at this time. The 
flow rate increased in the high-pressure mode to a value of 447 £/min (118 gal/min) at a 
time of 0,118 sec. At a time of 0,12 sec (see fig, 11(b)), the control initiated the transi- 
tion from the impact control limit force to the roll-out control limit force. During this 
transition, the control required a shift from the high-pressure mode of operation to the 
low-pressure mode at a time of approximately 0,139 sec. Fluid was then bled from the 
piston and reached a maximum flow rate of approximately -397 C/min (-105 gal/min). 

The flow rate decreased from this maximum until, at a time of 0,18 sec, the nose gear 
contacted the surface and the results were terminated. 

The control flow shown in this figure is the cumulative volume of fluid removed or 
injected into the piston at any instant of time. During the initial low-pressure mode of 
operation, approximately 0,156 fi (0,0055 ft^) of fluid was removed from the strut. This 
occurred at a time of 0.098 sec, where the transition from the low-pressure mode to the 
high-pressure mode of operation occurred. Fluid was then injected into the strut until the 
second operational transition occurred at a time of 0.139 sec. At this time all the fluid 
that had been removed from the strut had been returned and an additional amount of 
approximately 0.0255 i (0.0009 ft^) had been injected into the strut. Fluid was again re- 
moved from the strut, and at the time of nose -gear contact with the surface (0,18 sec), 
approximately 0,130 C (0,0046 ft^) had been removed. Consequently, after the impact 
phase, the fluid remaining in the gear was approximately 92 percent of the fluid in the gear 
prior to touchdown. 

Comparisons of computed results for the simply modified passive gear and the 
series -hydraulic control version of this gear are shown in figure 13 for a landing -impact 
simulation onto an extremely uneven, sinusoidally varying runway surface. Figure 13(a) 
presents wing force as a function of wing displacement. As in previous cases, the com- 
puted results show that the active control analysis is an energy-compatible procedure. 

Figure 13(b) shows time histories of shock-strut force and stroke for the landing 
impact onto the sinusoidally varying runway surface. The versatility of the active control 
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logic of the analysis may be seen by comparing the shock-strut -force results shown in 
this figure with those shown in figure 10(b) at the time of control activation. As discussed 
in connection with the vertical-drop simulation, control activation was initiated in the low- 
pressure phase, as anticipated during development of the control logic, and occurred dur- 
ing the fourth quarter of the first cycle of surface oscillation. In figure 13(b) control was 
initiated in the low-pressure phase during the first one -quarter cycle of surface oscilla- 
tion; however, the next one -half cycle of the surface oscillation was in phase with the wing 
motion and the shock-strut force for the passive gear decreased. As a consequence, the 
wing force decreased below the lower tolerance of the wing limit force and the high- 
pressure phase was required to maintain the wing force within the limits. Operation re- 
mained in the high-pressure phase until transition was required from the impact control 
limit force to roll-out control limit force. Since no experimental data are available on 
active control landing gears, this illustration of the versatility of the control logic in deal- 
ing with unexpected interactive effects increases confidence in the active control landing 
gear analysis. 

The maximum control flow rates encountered in this simulation were 432 C/min 
(114 gal/min) for the high-pressure phase of operation and -1154 f/min (-305 gal/min) 
for the low-pressure phase of operation. These large flow rates coupled with tire bounce 
(see fig. 13(b)) are undesirable. 

Time histories of vertical ground force and tire deflection for this simulation are 
presented in figure 13(c). The computed results for the active control gear show that 
upon control activation, ground force and tire deflection become greater than those for the 
passive gear. As previously noted in figure 13(a), the wing force for the passive gear is 
oscillatory because of the runway perturbations. It is desirable from the standpoint of 
fatigue not only to limit the magnitude of the forces applied to the upper body but to reduce 
oscillations as well. Consequently, with the control operating in the high-pressure mode 
the shock-strut force was increased to prevent the wing -force oscillation. As a result of 
this increased shock-strut force, the ground .force and tire deflection increased relative 
to those obtained for the passive gear. Presentations of active-gear results are discon- 
tinued at the time the tire leaves the surface. 

Landing roll-out simulations . - The effects of the series -hydraulic control on the 
damping characteristics of the modified oleo -pneumatic shock strut are shown in figure 14 
by comparing time histories of wing force and vertical ground force for the active and 
passive gears during a roll-out simulation. The simulations were made assuming that the 
mass of the airplane was supported by the main gears and the ground speed was 44.5 m/sec 
(146 ft/sec). 

Figure 14(a) presents computed time histories of wing force for the passive and 
active gears for the roll-out simulation. For the active gear, the wing force was con- 
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trolled about zero force with a tolerance of ±2134 N (±480 Ibf) and the wing force varied 
from zero to approximately -3114 N (-700 Ibf). At a time of approximately 0,08 sec the 
runway elevation dropped vertically to zero and the wing force for the passive gear in- 
creased very rapidly to a value of 4782 N (1075 Ibf), At a time of 2 sec, the wing force 
for the passive gear had damped to about 11 percent of the peak-to-peak value encountered 
during the first cycle. The active -gear simulation shows that the series -hydraulic control 
limited the wing force to the control limit force during the first two cycles, and at a time 
of 2 sec, the wing force had damped to about 5 percent of the peak-to-peak value encoun- 
tered during the first two cycles. 

Computed time histories of vertical ground force are presented in figure 14(b) for 
both the passive- and active -gear landing roll-out simulations. The peak-to-peak excur- 
sion of vertical ground force for the active gear, after the surface discontinuity, was 
greater than that for the passive gear. Subsequent peak-to-peak force oscillations, about 
the static ground force of 10,7 kN (2412.5 Ibf), were lower than those of the passive gear. 

On the basis of these comparisons of wing force and vertical ground force for both 
the passive and active gears, it appears that the damping characteristics of the series - 
hydraulic active control gear were equivalent to or better than those of the modified pas- 
sive gear for an excitation resulting from a vertically discontinuous, sinusoidal bump. In 
addition, the wing force encountered during traverse of the bump was reduced by 30 per- 
cent of the wing force experienced by the passive gear, 

CONCLUDING REMARKS 

A mathematical model of an active control landing gear has been developed and pro- 
gramed for operation on a digital computer. The principles upon which the model is based 
have also been discussed. The model has been evaluated for predicting airplane loads and 
motions encountered during simulations of a vertical drop and a landing impact of a rigid 
airframe with a landing gear incorporating a conventional (passive) oleo-pneumatic shock 
strut. Results have also been presented for the operation of a series -hydraulic active 
control gear for similar simulations of a vertical drop and landing impact and roll-out. 

The results obtained from this investigation are summarized in the following state- 
ments: The computed results from the subject mathematical model, simulating a vertical - 
drop test with a rigid airframe and a passive gear, are in good agreement with experimen- 
tal data and with results obtained from a multiple -degree -of -freedom takeoff and landing 
analysis (TOLA). Computed results for a landing -impact simulation, with a rigid airframe 
and passive gear, are in excellent agreement with results computed for the same landing - 
impact simulation employing TOLA, On the basis of these comparisons, it appears that 
the mathematical model of the active control landing gear is valid for the prediction of 
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airplane landing loads and motions for simulated symmetrical landings of a rigid airframe 
with passive main landing gears. 

Results of vertical-drop and landing -imp act simulations employing the mathematical 
model of the active control landing gear show that wing forces, relative to those encoun- 
tered with the modified passive gear, can be reduced substantially during the impact phase 
of the landing. The vertical -drop simulation for the active gear resulted in limiting the 
wing force to 75 percent of the force for the modified passive gear with a 15-percent in- 
crease in shock-strut stroke. The results of the landing-impact simulation for the active 
gear on a smooth runway indicated that the wing force could be limited to 80 percent of the 
force for the modified passive gear during the impact phase, with a 25 -percent increase in 
strut stroke relative to that encountered during passive-gear simulation. The results for 
a landing -impact simulation onto an extremely uneven runway also indicated that the wing 
force could be limited, during the impact phase, to 80 percent of the wing force encoun- 
tered with the modified passive gear. In this case, however, a 45-percent increase in 
shock -strut stroke relative to that encountered with the passive gear was required. These 
reductions in wing force could result in substantial increases in fatigue life of the 
structure. 

Comparison of computed results for the active and passive gears during a landing 
roll-out, in which the gears were forced into motion by a single, discontinuous, one- 
quarter cycle, sinusoidal bump, indicated that the series -hydraulic active control gear had 
damping characteristics equivalent to or slightly better than those inherent in the modified 
passive gear. For these simulations, the active gear reduced the wing force by 30 percent 
of the wing force generated by the passive gear. 

Control flow rates required by some of these simulations were very large and, if 
ultimately required for this type of active control gear, could necessitate advances in 
servo valve technology. However, it should be reemphasized that these results were for 
a simply modified version of a landing gear with a conventional oleo -pneumatic shock 
strut and, consequently, do not necessarily represent optimum results that could be ob- 
tained with a gear specifically designed to be compatible with the series-hydraulic control 
concept. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
December 18, 1975 
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APPENDIX A 


EQUATIONS OF MOTION 

This appendix presents the equations of motion for the physical system shown sche 
matically in figures 3 and 4 for the rigid-airframe option and the rigid -fuselage —elastic 
wing option. 


Rigid-Airframe Option 

The equations of motion for the rigid-airframe option in the body -axis system and 
the gravity-axis system are presented in this section. 

Airplane center -of -gravity (composite mass center) motion in body-axis system 


^c,b “ 


O.SLt sin a - 0.5(Df + E»t) cos a - 0.51>w + (0 + a) - Fq^^ sin (0 + a) 

- (mf + m^ + mh)g sin (0 + a) 


mf + m-v^7 + mh 


-0.5L^ - O.SLt cos a - 0.5(Df + D^) sin a - F^^gg + (mf + m^)g cos (0 + a) 


'c,b 


mf + m^y 


Airplane center -of -gravity motion in gravity-axis system 


Xc^g = Xq b cos (0 + a) + iZj, b sin (0 + a) 


zcjg = Zc,h cos (0 + O') - Xg^b sin (0 + a) 

Hub motion in body-axis system 

Xh,b = Xc^b + ^^2 + 0^di - 2 6s 


Fq z cos (0 + o) + Fq X sin (0 + a) + F^ gg + mbg cos (0 + o;) 

mb 
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Hub motion In gravity-axis system 

Xh,g = (0 + O') + z^^i, sin (0 + a) 

^h,g = ^h,b (0 + O') - xj^ t, sin (0 + a) 

Wing-mass -center motion in gravity-axis system 

^w,g ~ ^c,g ^ ^ cw '^cw + cos (^c\^ + 0; 

^Wjg ~ *^c,g ^ ^ cw cos (^cw "I" ~ sin 

Fuselage-mass-center motion in gravity-axis system 

Xf,g = '^C,g - sin (Xcf + 0) - 0^dcf cos (A-cf + 0) 

zf^g = ZQ^g - 0dcf cos (A_cf + 0) + 0^dgf sin (A-cf + 0) 

Wing-gear-interface motion in gravity-axis system 

^wG,g " ^c,g ^*^cwG (^cwG + ^) + ^ ^cwG (^cwG 

^wG,g ~ ^c,g ^^cwG (^cwG + ^) - ^^^cwG ("^cwG 

Shock -strut motion in body-axis system 

Sb= '^wG,g (0 + O') + x^G,g sin (0 + a) - zj^^g cos (0 + a) - sin (0 + a) 

Pitch angular acceleration about airplane center of gravity 

-0.53^dQ-j^ + 0.5D;fdQf + 0»5E)|^dQ^ + cos Of. - 0.5L^dQ-^^2 sin oi 

- 0.5D^dcw,z cos a - 0.5D^dcw,x sin a - Fa^ss^I + '^G,X^2 cos (0 + a) 

- ^G,Z^2 sin (0 + a) + Mn 
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Wheel angular acceleration about hub 

■■ ^ Ma + 
la 

Brake -control sensor -mass acceleration about hub 

is = 2Wn.)3 - ^ (0 - 4) + - /3) 

' '^/3,cr 

Rigid -Fuselage — Elastic -Wing Option 

The equations of motion for the rigid-fuselage — elastic -wing option in the body-axis 
system and the gravity -axis system are presented in this section. 

Airplane center -of -gravity (composite mass center) motion in body-axis system 

O.SLj sin a - 0.5(Df + Dt) cos a - 0.5Dw + cos (S + a) - Fq,Z sin (9 + o;) 

- (mf + m^ + mjj)g sin (9 + a) 

X . = -- — - 

mf + m^ + mji 


Zc,b - 


-0.5Lt cos Of - 0.5(Df + Dt) sin a - F^B + m^g cos (9 + a) 

mf 


- Fn cos 9 


Airplane center -of -gravity motion in gravity -axis system 
^c,g “ ^c,b (9 + a) + Zf, iy sin (9 + a) 

Z(.^g = Zg cos (0 + O') - X(j^ij sin (9 + a) 

Hub motion in body-axis system 

I^^G.X cos (0 + a) - Fq 2 sin (S + a) 

|- + m|^)g: sin (0 + ff) - 0 . 5 D^ 

|- O.SiDf + Dt) cos <y + 0.5Lt sin a 

Xh b = - — - + cos \eh + sin Xph 

n,u ^ ^ eii eii en on 

• * . Q 

+ sin \qq - 0 dee ^in A-ce 
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Prior to shock-strut stroking 

Fq^Z (6 + O') + Fq X sin (0 + a) 

+ + mh)g cos (0 + O') - 0.5 L^ + FwB 

2h,b = n, . + “deh sin \eh - cos X-eh 

Subsequent to shock-strut stroking 

^A,ss ^G,Z (0 + O') + Fq^x sin (0 + a) + mj^g cos (0 + or) 


Hub motion in gravity-axis system 

^h,g “ ^h,b (0 + a) + sin (0 + a) 

2h,g = ^h,b (0 + o) - sin (0 + a) 

Wing-mass -center motion in body-axis system 

Fq^X (0 + O') - Fq^ 2 sin (0 + a) 

- (mf + m^ + mh)g sin (0 + a) - 0.5F^B 

- 0.5(Df + D^;) cos a + 0.5L^ sin a 

.. L J . •* -o 

V, = — + Qfdp^ + Odcp sin - 0 dpp cos Xpp 

w,u ^ tjw oe oe oe 


Prior to shock-strut stroking 


z 


w,b 


^G, Z (0 + O') + Fq^^ sin (0 + O') + ^m^ + 
-0. 5L^ + ^wB 


m^ + mh 


+ ad0^ 


Subsequent to shock-strut stroking 

-^A,ss - + ^wB + + “) 

^w,b ~ ’ - + “dew 
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Wing-mass -center motion in gravity-axis system 
^W,g == ^w,b (0 + a) + 

z-^^g “ cos {0 + Q?) ~ x-^^b sin (6 + oi) 

Fuselage -mass-center motion in gravity-axis system 

%,g = ^c,g - ®dcf sin (\cf + 0) - 0^^^^ cos (A-cf + 

^f,g ~ ^c,g ■ ‘^°s (^cf + ®) + ^^dcf sin (A-cf + 0) 

Wing -gear -interface motion in body -axis system 

^wG,b ~ ^c,b ®^(*^ew + ‘^wwG) + ^"^cwG sin (^cwG - ®) + ^^^cwG ^°s (^cwG " “) 

-Fa,ss - 0-5Lw + FwB + mwg cos (0 + Q-) 

^wG,b ■- “ “ (*^ew + *^wwG) 

Shock -strut motion in body -axis system 
Sb = ^wG,b - ^h,b - “^^ss 

Motion of elastic axis at spanwise location of wing mass center: 

Prior to shock-strut stroking 

Fq 2 cos (0 + O') + Fq^x (0 + a) + (m^ + mh)g cos (0 + or) - 0.5L^ + F^g 

^ew.b ” " ' ““ — . — — 

’ + rtih 

Subsequent to shock-strut stroking 

-F^^gg - O.SL'^ + F^g + m-^g cos (0 + ot) 
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Motion of elastic axis at spanwise location of wing mass center relative to fuselage 

•2 

^er,b ^ew,b ” ^c,b ®^ce,b + ^ ^ce,b 

Rotational motion of wing about elastic axis and relative to fuselage: 

Prior to shock-strut stroking 

-Mgf + Fq^x^ss cos (0 + O') - Fg,Z^ss sin (0 + oi) - mhgZss sin (0 + a) 

+ Fa,z('^ew + ^wwG) cos (0 + O') + Fq^x(*^ 6W + ^wwO) sin (0 + a) 

^hS('^ew + '^wwG) cos (0 + oi) + m-v( 7 gdew + Ma 

CXq = ^ 

^ 2 2 

le + m^dew + la + mhdeh 

Subsequent to shock-strut stroking 

-Mef - FA,ss(dew + ^wwG) + m^gdew cos (0 + o;) + Fq,x^ss cos (0 + o;) 

- ^G,Z^ss sin (0 + O') - mhgfss sin (0 + a) + Ma 

= 

le + m^d^^ + la + mhd^j^ 


Wheel angular acceleration about hub 
" Ma + Mbr 


Brake -control -sensor acceleration about hub 


/3 = 2u)n, 13 


— ^(0 - /3) + cc.2 .(0 . /3) 
^/3,cr 


Pitch angular acceleration about airplane center of gravity 

-0.5Ltd3 + 0.5Dfd4 + 0.5Dtd4 + O.SL^ds cos a - O.SL^dg sin a 
-0.5E\vdcw,x cos a sin a - F^ sgdl + FQ,xd2 cos (0 + o') 
-FG,zd2 sin (0 + o;) + Mn 


39 



APPENDIX B 


DESCRIPTION OF LANDING GEAR 

The landing gear simulated in this paper was originally designed for a small air- 
plane having a gross mass of approximately 2268 kg (5000 Ibm). (See ref. 7.) The gear 
is a cantilevered type with a conventional oleo -pneumatic shock strut. The tire is a 
0.69-m (27 -in.) diameter type I (smooth contour) which was inflated to 221 kPa (32 psi). 
The mass of the landing gear is 68 kg (150 Ibm) and the unsprung mass is 59 kg (131 Ibm). 



The important geometric characteristics of the landing -gear shock strut are illus- 
trated in sketch (d), and the program input parameters are as follows: 


Shock strut: 

Pneumatic area, A 2 0.00535 m^ (0.05762 ft^) 

Hydraulic area, Aj 0.00437 m^ (0.04708 ft^) 

Primary orifice area, Aq 0.00005 m^ (0.00056 ft^) 

Pressurized pneumatic volume, V 2 0.0010 m^ (0.03545 ft^) 

Charging pressure, Pq a 299.9 kPa (6264 psfa) 

Bearing separation for fully extended shock strut, l\ 0.16828 m (0.5521 ft) 

Axial length from hub to lower bearing for fully extended 

shock strut, I 2 0.6785 m (2.22604 ft) 

Mass of airplane acting on each main gear, Wj 1094 kg (2411 Ibm) 

Unsprung gear mass, W 2 59 kg (131 Ibm) 
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Fully extended length of shock strut, ^ss,e 0.9674 m (3.174 ft) 

Friction coefficient (upper bearing), Mss 

Friction coefficient (lower bearing), Mss,Z 0.15 

Volume between piston and cylinder, Vo 0 (0 ft^) 

Area between piston and cylinder, A3 0 (0 ft^) 

Specific weight of hydraulic fluid, y^j 8226 N/m^ (52.36 Ibf/ft^) 

Dynamic viscosity of hydraulic fluid, /Xjj . . 0.00862 N-sec/m^ (0.00018 Ibf-sec/ft^) 

Mass density of hydraulic fluid, pjj 838 kg/m^ (1,626 slugs/ft^) 

Volume of hydraulic fluid in fully extended strut, Vjj .... -0.0014 (-0.05 ft^) 


Wheel and tire; 

Wheel flange diameter, df 

Width of wheel rim, Ww 

Unloaded diameter of tire, d 

Maximum width of undeflected tire, w 
Unloaded rated inflation tire pressure, Pj. 
Unloaded tire inflation pressure, Pq . . 

Pressure -rise parameter, k 

Vertical -force coefficient, 


. 0.39053 m (1.28125 ft) 
. . . 0.227 m (0.745 ft) 

0.686 m (nominal 2.25 ft) 
. . . 0.245 m (0.805 ft) 

482.63 kPa (10 080 psf) 
220.6 kPa (4608 psf) 

0.66 

0.02 
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Figure 1.- Variables considered significant for analyzing active control landing gear systems. 
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Figure 2.- Active control concepts for airplane landing gears. 
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Figure 5.- Operational phases of series -hydraulic active control gear. 
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O Experiment (ref. 7) 
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(a) Time histories of shock-strut force and stroke. 


Figure 7.- Comparison of computed results obtained from ACOLAG and TOLA with experimental data for a vertical 
drop onto a flat stationary surface. Q = 0^; = 2.7 m/sec (8.8 ft/sec); Xc^g = 0 m/sec (0 ft/sec). 
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(a) Wing force as a function of wing displacement. 

Figure 9.- Comparison of computed results for modified passive and active gears for a vertical drop 
onto a flat stationary surface. 6 = 0°; Zc,g = 2.7 m/sec (8.8 ft/sec); Xc,g = 0 m/sec (0 ft/sec). 
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(c) Time histories of vertical ground force and tire deflection. 
Figure 9.- Concluded. 




(a) Wing force as a function of wing displacement. 


Figure 10.- Comparison of computed results for modified passive and active gears for a vertical drop onto 
a vertically oscillating surface. 9 = 0^; = 2.7 m/sec (8.8 ft/sec); xc,g = 0 m/sec (0 ft/sec). 
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(a) Wing force as a function of wing displacement. 

Figure 11.- Comparison of computed results for modified passive and active gears for a landing impact onto 
a smooth runway surface. 0 = 10.5^; Z(* g = 2.7 m/sec (8.8 ft/sec); xc^g = 45.7 m/sec (150 ft/sec). 
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Figure 11.- Continued. 
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Figure 12.- Time histories of control flow rate and control flow for landing impact of active gear onto a smooth 
runway surface. 9 = 10.5^; 2.7 m/sec (8.8 ft/sec); = 45.7 m/sec (150 ft/sec). 
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(a) Wing force as a function of wing displacement. 

Figure 13. - Comparison of computed results for modified passive and active gears for a landing impact onto 
sinusoidal runway surface. 6 = 10.5°; = 2.7 m/sec (8.8 ft/sec); Xc^g = 45.7 m/sec (150 ft/sec). 
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(b) Time histories of shock-strut force and stroke. 
Figure 13.- Continued. 
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(a) Time histories of wing force. 


Figure 14.- Comparison of computed results for active and passive gears for a landing roll-out after traversing 
vertically discontinuous bump. 0 = OO; Zc g = 0 m/sec (0 ft/sec); Xc g = 44.5 m/sec (146 ft/sec). 
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